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Results of modeling the spectra of two supernovae SN 2008D and SN 2006aj related to the 
X-ray flash XRF 080109 and gamma-ray burst GRB / XRF 060218, respectively are studied. 
The spectra were obtained with the 6-meter BTA telescope of the Special Astrophysical 
Observatory of the Russian Academy of Sciences in 6.48 and 27.61 days after the explosion 
of SN2008D, and in 2.55 and 3.55 days after the explosion of SN2006aj. The spectra were 
interpreted in the Sobolev approximation with the SYN0W code. An assumption about the 
presence of envelopes around the progenitor stars is confirmed by an agreement between the 
velocities of lines interpreted as hydrogen and helium, and the empiric power-law velocity 
drop with time for the envelopes of classic core-collapse supernovae. Detection of a P Cyg 
profile of the H/3 line in the spectra of optical afterglows of GRBs can be a determinative 
argument in favor of this hypothesis. 

1. INTRODUCTION explosion energy. If the helium envelope is lost 

and the spectra show only the traces of heavier 

The modern classification of core-collapse su- elements (O, Mg, Si, S, Ca, Fe), then the super- 

pernovae includes several subtypes. If the spec- nova is classified int ° the Ic type. Different sub- 

tra demonstrate hydrogen lines, then the super- ^ es of supernovae vary by different explosion 

nova belongs to the II type. Its light curve ener g ies > and > consequently, by the velocities of 

shape is determined by the hydrogen recom- matter ex P ansion in * he envelopes, 

bination features in the envelopes. If a star Since in some cases the classification depends 

loses its hydrogen envelope, but the helium shell on the phase of the supernova at the moment 

is still present in the spectrum, then the star of observation, some intermediate types are in- 

is classified into the lb type. The impacts of traduced: Ib/c (appearance of weak helium lines 

hydrogen recombination are absent in its light in late spectra), lib (disappearance of hydrogen 

curve, and its shape is determined by the heat- lines in late spectra). The classification of core- 

ing of the emitting matter as a result of the collapse supernovae can be conceived as a mani- 

56 Ni — > 56 Co — > 56 Fe decay, depending on the fold of intermediate subtypes, depending contin- 
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uously on the initial parameters of the progenitor 
star. 

This paper discusses the manifestations of 
stellar-wind envelopes around the progenitor 
stars of an X-ray flash XRF 080109 and a 
gamma-ray burst with a strong X-ray compo- 
nent GRB / XRF 060218. Some signs of super- 
novae (SNe) were noticed in the spectra and light 
curves of both events, which allowed studying 
these phenomena from the very start of the ex- 
plosion (unlike most SNe that are detected at the 
moments close to the "classical" maximum of the 
light curve). The study of this kind of phenom- 
ena makes it possible to approach the solution 
to the supernova explosion mechanism and the 
origin of gamma-ray bursts. 

On January 9.57, 2009 (UT) the X-ray tele- 
scope XRT on board of the Swift space observa- 
tory registered an X-ray flash XRF 080109. The 
X-ray component was detectable for about 15 
minutes [1,2]. The Burst Alert Telescope (BAT) 
did not detect any gamma-ray quanta. The pas- 
sage of a shock wave through the envelope was 
observed in multi-color light curves during sev- 
eral days [2, 3]. The afterglow light curves and 
spectra distinctly showed signs of the supernova, 
designated as SN2008D. 

Almost two years earlier (on February 18.149, 
2006) the same space platform Swift detected 
a peculiar gamma-ray burst with the SN signs 
in the afterglow spectra and light curves of 
GRB /XRF 060218/ SN2006aj. It is the first 
phenomenon for which we observed the passage 



of a shock wave through the wind envelope sur- 
rounding a massive SN progenitor star, which 
manifested itself as a thermal component in the 
X-ray spectrum observed during the first 2 hours, 
and then powerful UV burst with the maxi- 
mum observed in 11 hours after the burst [4, 5]. 
The effect of the shock wave passage is seen in 
the optical light curve as a relatively short peak 
several days of duration [5, 6]. 

In both cases the moment of the X-ray flash 
and/or gamma-ray burst registration was taken 
as the beginning of the SN explosion. 

The structure of the paper is as follows: spec- 
tral observations and preliminary data process- 
ing are described in section 2; the interpretation 
of the obtained spectra with the SYNOW code is 
laid out in 3; a comparison of the obtained mea- 
surements of the photosphere and envelope ex- 
pansion velocities with the analogous velocities of 
SNe that did not show any relation with gamma- 
ray bursts or X-ray flashes is in section 4; our 
conclusions are discussed in section 5. 

2. SPECTRAL OBSERVATIONS 

Optical afterglows of XRF 080109 and GRB / 
XRF 060218 were observed with the SCORPIO 
Optical Reducer mounted in the primary focus of 
the 6-m BTA telescope of the Special Astrophys- 
ical Observatory of the Russian Academy of Sci- 
ences (SAO RAS). The VPHG550G grism cov- 
ering the spectral range of 3500-7500 A with a 
resolution (FWHM) of 10 A was used as a dis- 
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persing element. 

We obtained two spectra of SN 2008D: on Jan- 
uary 16 and February 6, 2008 (6.48 and 27.61 
days after the explosion) . For SN 2006aj two 
spectra were obtained as well: on February 20 
and 21, 2006 (2.55 and 3.55 days after the ex- 
plosion). The data processing was standard and 
included: subtraction of the electronic zero (an 
additive component of the total signal produced 
by the CCD chip), flat-field correction (com- 
pensation of the CCD chip sensitivity irregular- 
ities), wavelength calibration with the aid of a 
comparison spectrum of a Ne-Ar lamp, atmo- 
sphere extinction correction, and calibration by 
the absolute flux with the use of observations of 
a spectrophotometric standard performed every 
night. Apart from that, the spectra of SN2008D 
were subtracted for the contribution of the host 
galaxy, the spectral distribution of which is con- 
structed from the regions located in immediate 
proximity to the supernova. 

Then the observed spectra were corrected 
for the extinction in the Galaxy according 
to the dust distribution maps [7]. When 
taking the extinction into account, the dust- 
screen model was accepted, where the ex- 
pression to account for extinction is of the 
form F int {\) = F obs {\)W 0Axk ^ E( - B - v \ where 
FintW and F b s (\) are the emitted (without ab- 
sorption) and observed fluxes, respectively. The 
Milky Way extinction curve k(\) was taken from 
the paper [8]. The spectra of SN2008D were 
as well corrected for the extinction in the host 



galaxy according to the data from [3]. The ab- 
sorption in the host galaxy of SN 2006aj is neg- 
ligible [9]. 

Before the interpretation, the spectra were 
transferred to the reference frames, associ- 
ated with the gamma-ray burst or the X-ray 
flash (z = 0.0331 for GRB / XRF 060218 and 
z = 0.007 for XRF 080109, see Fig. 1). The red- 
shifts were estimated from the shift of galactic 
emission lines, and are consistent with data pub- 
lished in the literature. 

3. COMPARISON OF THE OBSERVED 
AND SYNTHETIC SPECTRA 

To interpret the observed spectra in detail, 
we applied the multi-parametric SYNOW code [10], 
which was previously used for the analysis of 
spectra of core-collapse supernovae [11-13]. The 
code algorithm is based on the following assump- 
tions: spherical symmetry; homologous expan- 
sion of layers (v ~ r); sharp border of the photo- 
sphere emitting a black-body spectrum and as- 
sociated at early stages with a shock wave. 

The code is used to identify the lines and find 
the expansion velocities of the layers in which 
they are formed. The P Cyg line profiles ob- 
served in the spectra of supernovae and modeled 
with the code are then divided into two types 
according to their shape: 

• the lines formed in the layer undetached 
from the expanding photosphere; 
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Figure 1. Spectra of the supernova SN2008D obtained with the BTA on January 16 (black line) and 
February 6 (grey line), 2008 in 6.48 and 27.62 days after the explosion, respectively. 



the lines appearing in the layer detached 
from the photosphere. 



The version of undet ached layers was used for 
the first spectrum of SN2006aj. A combined 
version in which lighter ions are detached and 
the heavy ions are not detached from the pho- 
tosphere was applied for the second spectrum of 
SN 2006aj and for both spectra of SN 2008D. The 
version choice was determined by fitting the pa- 
rameters of every ion to the observed spectral 
features. 

The model parameters of the input file and 
examples of interpretation of the SNe spectra are 
described in detail in the papers by code devel- 
opers [10-13]. 



3. 1 . Modeling of SN 2008D Spectra 

The main absorption features are noticeable 
in both spectra (see Figs. 2 and 3) as P Cyg 
profiles of He I, Fell, OI and, presumably, HI 
(a spectral feature near 6200 A) lines. The main 
models containing HI as a candidate describing 
the above-mentioned feature are shown by the 
thick black line, and the observed spectra— by 
the thin grey noisy line. 

The first spectrum is best described 
by the model with the following pa- 
rameters: photosphere temperature 
T55 = 8700 K, velocity of the photosphere 
Vphot = 17000 km/s, maximum matter expansion 
velocity in the envelopes V max = 70000 km/s. 
The parameters of the best model for the 
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Figure 2. The spectrum of the supernova SN2008D obtained with the BTA on January 16, 2008 in 6.48 
days after the explosion (in the reference frame linked with the object, z = 0) is shown by the thin grey 
noisy line; the thick black line marks the model spectrum which contains HI; the model spectrum containing 
C II instead of hydrogen is marked by the thin dashed line; the model spectrum with Si II instead of 
hydrogen is marked by the solid thin black line. 



second spectrum are as follows: = 6000 K, 
Vphot = 8500 km/s, V max = 70000 km/s. 

The temperature of the photosphere was de- 
termined by the slope of the observed spectrum. 
Since the shape of the spectra is very differ- 
ent from the Planck curve, the ions, describing 
the main spectral features were immediately in- 
cluded in the calculation. The photosphere ve- 
locity was adopted as equal to the velocity of the 
ionized iron lines [13]. The maximal velocity at 
the upper boundary of all envelopes V max was 
estimated from the blue part of the P Cyg line 
profiles. Optical depth of the reference line [13] 
was selected in a way to best describe the ab- 



sorption components of the P Cyg profile. The 
SYN0W code allows applying one of the two laws of 
the matter density decrease with distance: power 
(r(r) oc v(r)~ n , where v ~ r), and exponential 
(r(r) oc exp(— v (r)/V e ), where V e is the charac- 
teristic velocity of the layer, which determines 
the size of the envelope and the width of the ab- 
sorptions in the observed spectrum). We used 
the exponential law in our calculations. The pa- 
rameters of the main models that take hydrogen 
into account are presented in Table 1. 

As an alternative interpretation of the 
spectral feature near 6200 A, we considered the 
model calculations of synthetic spectra contain- 



6 



i — r~r~i — i — 1 1 1 i i — | — i — r~r~i — i — i — r~r~i — | — i — r~r~i — i — i — r~r~i — | — i — r~r~i — i — i — i — r~i — | — i — i — r~i — i — i — i — r~r 




Figure 3. The spectrum of the supernova SN2008D obtained with the BTA on February 6, 2008 in 27.62 
days after the explosion. The markings are the same as in Fig. 2. 

Table 1. Parameters of the models which best describe both spectra of SN2008D. The parameters correspond to 
the models, marked with thick black lines in Figs. 2 and 3; r is the optical depth of the reference line of every ion; 
Vmin is the minimal layer velocity of every ion; V e is the characteristic velocity in the applied relation between the 

optical depth and the expansion velocity r(r) ~ exp(— v(r)/V e ), where v ~ r; V p hot is the photosphere velocity; 
Vmax is the velocity at the upper boundary of all envelopes; Ttb is the black body temperature of the photosphere, 

forming a continuous spectrum. 



Parameters 


First spectrum (6.48 days) 


Second spectrum (27.62 days) 


Ions 


HI 


He I 


Fell 


01 


HI 


He I 


Fell 


01 


T 


0.4 


0.5 


1.0 


0.5 


0.5 


3.0 


2.0 


0.3 


V min , km/s 


23000 


17000 


17000 


17000 


15000 


10500 


8500 


10500 


V e , km/s 


8000 


4000 


3000 


4000 


3000 


1000 


3000 


4000 


Vphot, km/s 


17000 


8500 


V m ax, km/s 


70000 


70000 


T bb ,K 


8700 


6000 



ing the lines of Sill and CII instead of the HI supernovae spectra [13, 14]. In Figs. 2 and 3 the 
lines. These ions were frequently mentioned thin solid lines correspond to the models with 
in the literature on modeling the type Ib-c silicon taken instead of hydrogen: r = 0.005, 
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V min = 17000 km/s, V e = 10000 km/s for 
the first spectrum (all the other parame- 
ters are identical to those in the model with 
hydrogen), and r = 0.5, V m i n = 8500 km/s, 
V e = 5000 km/s for the second spectrum. 
The models containing the Sill lines 
(Kest = 6347 A) instead of HI lines are 
velocity-limited by the photosphere and can 
not describe the observed spectral features near 
6200 A according to the modern concepts on the 
radial stratification of elements [13]. 

The absorption feature near 6200 A can 
as well be explained by the presence of ion- 
ized carbon CII with the following model 
parameters: r = 0.005, V m in = 24000 km/s, 
V e = 10000 km/s for the first spectrum, 
and r = 0.0008, V min = 16000 km/s, 

V e = 3000 km/s for the second spectrum 
(the thin dashed lines in Figs. 2 and 3). 

Thus, Sill can be excluded from the calcula- 
tions of these two spectra of SN2008D. CII line 
(Kest = 6580 A) remains a possible alternative 
for the description of the absorption feature near 
6200 A. The above mentioned feature can possi- 
bly be formed by a blend of the Ha and C II lines. 
It is obvious that the final conclusion requires ob- 
servations with a higher spectral resolution. 

3.2. Modeling of SN2006Aj Spectra 

The spectra of SN2006aj obtained with the 
BTA in 2.55 and 3.55 days after the gamma- 
ray burst were modeled with the SYN0W code. 



The first observed spectrum is best modeled 
with the photosphere expansion velocity of 33000 
km/s. This parameter is within measure- 
ment errors of the photosphere expansion ve- 
locity, measured from the Swift/XRT/UVOT 
data ((2.7±0.8)xl0 4 km/s [4]). At the mo- 
ment we started our spectral observations 
(roughly straight after the last observation of 
the UV flash with the Swift/UVOT), the photo- 
sphere velocity remained within approximately 
the same limits. Moreover, the early spec- 
trum shows a broad continuum depression at 
5900-6300 A and an almost imperceptible ex- 
cess at 6300-6900 A, which is best described by 
the Ha P Cyg profile with the same velocity of 
33000 km/s. For the fitting we used different 
models of the case where the layers are unde- 
tached from the photosphere, but the velocities 
of the photosphere and ions remained equal to 
33000 km/s. The model parameters, the fitting 
of the observed spectrum, and detailed comments 
are given in [6]. 

The velocities best describing the second 
spectrum (3.55 days after the burst) are within 
the limits 18000 km/s < {V phot , V min } 24000 
km/s. The parameters of the model spectrum de- 
scribing the observed spectrum, and the selected 
values of model parameters are laid out in [6] . 

Using the best models of the SN2006aj and 
SN2008D spectra, we measured the hydrogen 
and helium envelope velocities, and the veloci- 
ties of photosphere expansion for the respective 
observation epochs (see Table 2). The earliest 
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Figure 4. Comparison of light curves of the following supernovae in the V-band: SN 1983N [15-17]; 
SN1984L [17]; SN 1999dn [18]; SN2000H [19]; SN2008D [3]; SN2006aj [6[. The thick grey line shows the 
light curve of SN 2006aj without scaling along the time axis; the thick black line shows the same light curve 
scaled along the time axis with a factor of 1.7. The SN2008D light curve is shown with the black line 

connecting the points with bars. 



spectra of SN 2006aj demonstrate one the high- 
est observed expansion velocities of the order of 
0.1 of the speed of light. The teams of Kunugasa 
and Sahu ([21, 22]) describe similar expansion ve- 
locities observed in the early spectra of type Ic 
supernovae SN 2002ap and SN 2007ru, which in- 
dicates the presence of not just a unique object 
but rather a class of events capable of showing 
such high velocities. 

We compared the velocities obtained for 
SN2008D and SN20006aj with an empiric law 
of photosphere and envelope velocity decrease for 
11 type lb supernovae from [11] using the scal- 



ing factor for SN 2006aj that was obtained based 
on a comparison of the light curves of the super- 
novae sample. This comparison is presented in 
Figs. 5 and 6. 

4. PHOTOSPHERE AND ENVELOPE 
VELOCITIES 

Light curves of different type lb, Ic and Ib- 
c supernovae mainly differ in absolute fluxes, 
widths of bell-shaped peaks near the main max- 
ima, and their behavior at late phases of the 
flash. However, the light curve shapes of dif- 
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Table 2. Velocity variation of the hydrogen and helium envelopes and photosphere expansion of the 
supernovae SN 2006aj and SN 2008D after the explosion. The values in brackets indicate time after the 
maximum. The time of the maximum for SN 2006aj was determined as 10.4 days after the explosion (based 
on the V-band light curve) [9], and 19 days for SN2008D [20]. 
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Figure 5. Variation of photosphere velocities, measured from the Fe II lines [11]. Photosphere velocities 
measured from our spectra are denoted as triangles for SN2008D and squares for SN2006aj. The scaling 
factor 1.7 was applied to the SN2006aj photosphere velocity. 



ferent supernovae near the maximum are similar SN 2008D and SN 2006aj in different bands with 

to each other. Valenti et al. [14] compared the the light curves of SNe from [11]. In the case 

light curves of SNe of different types via seal- of SN2006aj, we applied the scaling factor 1.7 

ing along the time axis. We applied an identi- (see Fig. 4). A comparison of light curves of the 

cal method when comparing the light curves of SNe not observed in the V-band with the data 
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Figure 6. Minimal velocities from He I lines (filled squares for the case of layers undetached from the 
photosphere, filled diamonds for the detached lines), and HI lines (empty circles, always in detached layers) 
[11]. The curve represents the photosphere velocity drop power law from Fig. 5. Minimal velocities of He I 
lines for SN2008D are indicated by triangles up, HI line velocities — by triangles down. The velocities of 
hydrogen and helium lines for SN2006aj are indicated by squares. The scaling factor 1.7 was applied to 

SN 2006aj envelope velocities. 



available in the B- and R-bands confirmed the 
accuracy of our choice of the scaling factor for 
the SN 2006aj light curve. 

5. CONCLUSIONS 

An agreement between the velocities of 
photosphere expansion and those of the he- 
lium/hydrogen envelopes of the SNe associ- 
ated with the X-ray flash XRF 080109 and the 
gamma-ray burst GBR 060218 argues for similar 
explosion dynamics of gamma-ray bursts and X- 



ray flashes. 

A comparison with analogous velocities of 
classical type lb supernovae not accompanied by 
gamma- or X-ray bursts may indicate common 
properties of the progenitors of these two classes 
of phenomena, and similar explosion mecha- 
nisms. 

Previously, an idea about the presence of hy- 
drogen in the spectra of type lb and Ic super- 
novae was suggested [13, 23]. The idea was con- 
firmed in our paper [6] dedicated to the study 
of early spectra of SN2006aj. 
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Branch et al. [23] list the main candidates 
for the description of a spectral feature near 
6200 A observed in the early spectra of super- 
novae. Apart from the hydrogen line Ha A6563 
these are: Sill A6355 spaced from Ha at 9500 
km/s, Nel A6402 at 7360 km/s, and only CII 
A6580 is spaced from the hydrogen line at 777 
km/s, which can be within velocity measurement 
errors. Candidate velocities must agree with the 
photosphere velocity determined from the heavi- 
est ions (e.g., from the ionized iron). Reasoning 
from the results of modeling the SN 2006aj and 
SN2008D spectra, it can be concluded that the 
most possible sources of absorption near A6200 
are Ha A6563 and CII A6580. To choose be- 
tween these candidates we compared the velocity 
variability of hydrogen envelopes of the two SNe 
under study with the classical supernovae. The 
results of comparison of the velocity evolutions 
of different supernovae argue for the existence 
of hydrogen-containing envelopes around massive 
progenitor stars that appear before the explosion 



due to more or less powerful stellar winds. 

Detection of P Cyg profile absorption com- 
ponents of other Balmer lines, in particular H/3 
A4861, can become a decisive argument in favor 
of hydrogen presence in the envelope. Since the 
optical depth of H/3 line is smaller than that of 
Ha, there arises a problem of its detection in a 
noisy spectrum among numerous lines of heavy 
metals, particularly, Fell. Detailed long-term ob- 
servations of the evolution of subtle details in the 
spectra obtained with a high signal/noise ratio 
can solve this problem. 
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